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Abstract:  Presented  in  this  memorandum  are  the  theoretical  modifications 

to  the  radial  equilibrium  equation  and  the  computerized 
adaptations  required  to  approximate  the  effects  of  skewed 
blading  in  the  current  working  version  of  the  streamline 
curvature  method  computer  program  (SCM-IV) . The  procedure 
is  illustrated  by  the  analysis  of  a pumpjet  that  is  fitted 
with  a skewed  rotor.  The  results  of  this  analysis  which 
utilizes  the  data  generated  by  PRESKEW  show  an  improved 
agreement  with  experimental  data  relative  to  SCM-IV 
calculations  prior  to  the  skew  modification. 
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Introduction 

The  initial  version  of  the  streamline  curvature  method  computer  program 
(SCM-I,  Ref.  [a])  for  the  analysis  of  axisymmetric  turbomachine  flows  was 
applicable  only  tc  flow  regions  outside  of  blade  rows.  To  design  skewed 
blading,  a method  was  needed  to  approximate  the  flow  effects  within  the 
blade  rows.  Bruce,  in  Appendix  B,  has  shown  that  the  radial  equilibrium 
equation  for  application  within  blade  rows  developed  by  Smith  in  Ref.  [b] 
was  amenable  to  solution  by  streamline  curvature  techniques.  McBride  in 
Appendix  C,  has  offered  additional  physical  supporting  arguments  relevant 
to  the  effects  of  blade  skew. 

This  paper  will  summarize  the  work  of  the  three  previously  mentioned 
authors  and  present  the  approach  used  to  successfully  modify  SCM-I.  (The 
current,  working  version  of  the  streamline  curvature  method  program,  SCM-IV, 
to  be  documented  at  a later  date,  includes  approximations  of  shroud  inlet 
losses,  blade  secondary  flow  losses,  and  the  effects  of  preswirl,  counter- 
rotation, and  blade  skew.) 

PRESKEW  is  an  independent  computer  program  which  develops  the  required 
data  for  the  utilization  of  the  "skew  option”  in  SCM-IV.  At  this  time,  it 
should  be  emphasized  that  PRESKEW  and  SCM-IV  represent  a first  approximation 
of  the  blade  skew  effects.  However,  the  comparative  results  presented 
later  in  this  paper  demonstrate  encouraging  agreement  between  computed  and 
experimental  data. 

Analysis 

When  heat  addition  and  compressibility  effects  are  neglected,  the 
radial  equilibrium  equation  for  application  within  blade  rows  has  the 
following  form: 


V2  V2  3V  V 

13P  0.  i m . . EL  . z ► . DA  v 

~ -r—  = — + cos  4>  n V sin  4>  -s — + -r-  tan  q>  — + F 

p9rr  R m dm  A Dz  r 

m 


where 


p *=  fluid  mass  density,  slugs/ft3 
P = local  fluid  static  pressure,  psfa 
r = radial  coordinate,  ft 
z = axial  coordinate,  ft 
m = meridional  coordinate,  ft 


4>  = meridional  flow  angle  with  respect  to  the  centerline,  radians 
VQ  “ peripheral  velocity  component,  ft/sec 
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Vffi  = meridional  velocity  component,  ft/sec 

V = axial  velocity  component,  ft/sec 
z 

R = radius  of  curvature  of  the  meridional  streamlines,  ft 

m 

X = effective  area  coefficient 

Fr  *=  radial  component  of  the  body  force  per  unit  mass,  lbs/slug 
D = differential  change  following  particle. 


Equation  (1)  differs  from  the  radial  equilibrium  equation  that  is 
applied  outside  of  blade  rows  by  the  presence  of  the  last  two  terms  on 
the  right-hand  side  of  the  equal  sign.  The  first  of  these  two  terms  is 
a function  of  the  blade  blockage.  While  the  Fr  term  depends,  inpart, 
on  the  blade  skew.  With  reference  to  Figure  1,  the  local  blade  skew 
angle,  £,  is  defined  as  the  angle  formed  in  the  R-0  radial  plane  between 
the  tangent  line  at  a point  to  the  blade  mean-thickness  line  and  a 
radial  line  through  the  same  point. 

The  effective  area  coefficient,  X,  is  given  by 


X 


s 


V 


(2) 


where 


N = the  number  of  blades 

0S  * the  peripheral  angular  coordinate  of  the  blade  suction 
surface,  radians 

0p  *=  the  peripheral  angular  coordinate  of  the  adjacent  blade 
pressure  surface,  radians. 


Although  Bruce  has  indicated  the  significance  of  the  DA/Dz  term  in 
Appendix  B,  most  ARL  blades  are  thin  and  parallel  at  a given  section. 


and 


DA 

Dz 


is  neglected 
1.0 

(3) 

0.0 

(4) 
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Therefore,  for  the  PRESKEW  approximation  the  first  of  the  additional  terms 
is  omitted. 

Bruce  and  McBride  have  both  shown  that  the  radial  component  of  the  body 
force,  Fr  , can  be  expressed  as 


AP  tan  £ 


r pr(0s  - 6p) 


(5) 


where  AP  is  the  local  static  pressure  difference  across  the  blade  pressure 
and  suction  surfaces. 

If  a local  static  pressure  coefficient  across  the  blade  is  defined  as 


C'  = 


AP 


P 1/ 2pV  2 


where  V is  a free  stream  or  reference  velocity  in  ft/sec,  then  the  dimension- 
00 

less  radial  component  of  the  body  force  per  unit  mass,  Fy,  can  be  written  as 


N tan  £ C* 

F — E- 

2 | r 4ttX  (r/R  ) 


(6) 


where  RQ  is  a reference  radial  distance  such  as  the  maximum  body  radius  or 
propeller  tip  radius.  When  the  radial  derivative  of  Bernoulli’s  Equation 


13P  l9Pt  13V’ 


13VI 


JL  w*  J-  t -L  m x 

p 3r  p 3r  2 3r  2 3r 


(7) 


where  Pt  is  the  fluid  total  pressure,  is  combined  with  Equations  (1)  and  (6) 
in  a dimensionless  manner,  the  "skewed"  version  of  the  radial  equilibrium 
equation  for  computations  within  blade  rows  results: 


3 CM 
3Y 


Cos<j>  Sin<ft  3 CM 
RC  “ CM  M 


CM  = — 
CM 


1 3CPT  CU  3 (Y»CU)  _ 

2 3Y  ” Y 3Y  Y 


(8) 


where 


Y 

M 


r/Rc 

m/R 


dimensionless  radial  distance 
dimensionless  meridional  distance 
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RC 

- R /R 

m o 

CM 

= V /V 

m oo 

CU 

- VVo° 

P - p 

CPT 

t 

1/2PVJ 

dimensionless  radius  of  curvature 
meridional  velocity  ratio 
peripheral  velocity  ratio 

total  pressure  coefficient 


\ 


i 

I 


To  apply  Equation  (8)  within  the  blade  row  the  changes  in  the  angular 
momentum  coefficient,  AMq  , and  total  pressure  coefficient,  ACPT  , from  the 
blade  leading  edge  to  the  computing  station  must  be  specified  as  a function 
of  Y. 


AM0  = Y2  • CU2  - Yx  • CUX 


(9) 


where  the  subscripts  1 and  2 refer  to  the  blade  leading  edge  and  the  computing 
station,  respectively.  From  Euler's  turbomachine  equation 


ACPT 


(u/vj. 


cu2  - (u/vj1 


(10) 


where  U/V  is  the  dimensionless  rotor  wheel  speed  and  is  defined  by 
00 


u/vto  = 


ttY 

J 


(ID 


V 

oo 

J = — = the  advance  ratio 
nD 


(12) 


where 

n = the  rotational  speed  of  the  rotor,  rps 

D = 2R  = a reference  diameter,  ft. 
o 


When  J has  been  specified  AMq  and  ACPT  are  functions  only  of  CU  for  a 
particular  value  of  Y at  the  computing  station.  For  blades  designed  by  the 
mean  streamline  method  of  Wislicenus,  Reference  [c],  AVg/V0t  and  AP/AP  are 
empirically  specified  as  a function  of  the  axial  blade  length  ratio  x/ £, 
where  at  any  blade  cylindrical  section 

x “ the  axial  distance  from  the  blade's  leading  edge,  ft 


_ A 
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the  axial  projection  of  the  blade  chord  at  the  given  section, 
ft 


the  change  in  peripheral  velocity  from  the  blade's  leading  edge 
to  the  computing  station,  ft/sec 

the  total  change  in  peripheral  velocity  from  the  blade's  leading 
edge  to  trailing  edge,  ft/sec 

the  local  static  pressure  difference  across  the  blade  at  the 
computing  station,  psfa 

the  mean  static  pressure  difference  for  the  entire  blade  section, 
psfa. 


Once  these  data  have  been  specified.  Equations  (9)  and  (10)  can  be  evaluated 
from  the  known  CU^  versus  Y distribution  by  employing  the  following 
relationship: 


cu2  = 


The  AP/AP  distribution  is  used  to  calculate  the  radial  body  force  para- 
meter, Fy.  If  £ is  specified  as  a function  of  Y,  Fy,  Equation  (6),  for  any 
cylindrical  section  becomes  a function  of  x/Jl.  For  the  fluid  element  shown 
in  Figure  2 which  has  dimensions  t by  Z by  dr,  the  peripheral  component  of 
the  pressure  force,  dF0  , must  be  equal  to  the  change  of  momentum  in  the 
peripheral  direction: 


dF0  = AP£dr  = pVgtdr  AV0t 


where  Va  is  an  average  meridional  velocity  for  the  particular  blade  section. 
If  Equation  (14)  is  rearranged  and  divided  by  l/2pV 2,  C'  can  be  calculated 
from 


— /V 

AP  = 2 _a 

l/2pV  2 


l/2PVa 


9 March  1977 
ALT : jep 


r 


where 


t_  = 2ffr 
£ = N£ 


(17) 


Thus,  by  combining  Equations  (6),  (9),  (10),  (11),  (12),  (13),  (15), 
and  (16)  , the  necessary  parameters  for  the  utilization  of  Equation  (8) 
within  the  blade  row  can  be  calculated.  PRESKEW  was  developed  to  compute 
Fy,  AMQ,  and  ACPT  as  a function  of  Y at  the  computing  station. 


Computerization 


[The  logical  operation  of  PRESKEW  is  documented  by  the  functional 
diagram  of  Appendix  A.  The  following  discussion  will  summarize  the 
operation  of  the  program  and  explain  the  approximations  and  procedures 
used  to  calculate  Fy,  AM0 , and  ACPT  as  a function  of  Y for  use  in  SCM-1V. 
Most  of  the  data  specification,  interpolation,  and  integration  is  performed 
by  the  spline  curve  procedures  documented  by  Davis  in  Reference  [d]. 


To  avoid  the  convergence  problems  inherent  to  streamline  curvature 
techniques  when  using  high  aspect  ratio  computing  stations,  one*  axial 
station  within  the  blade  row  was  chosen  for  the  radial  specification  of  the 
skew  effects.  The  blade  planform  shape  is  represented  by  four  inter- 
secting straight  lines  as  depicted  in  Figure  3.  The  blade  skew  angle, 

£,  is  assumed  to  be  constant  for  the  current  version  of  PRESKEW.  This 
assumes  that  the  blade  mean-thickness  line  is  described  by  a logarithmic 
spiral  as  illustrated  in  Figure  1. 

Once  the_input  data  have  been  read,  computation  begins  by  approx- 
imating the  (Va/V  ) distribution  as  a function  of  Y along  a radial  line 
midway  between  the  two  computing  stations.  This  calculation  is  based 
on  the  continuity  principle  and  integration  of  the  reference  meridional 
velocity  profile. 

As  the  computing  station,  the  radial  computation  points  (values  of  Y) 
and  the  projected  local  blade  section  length  are  computed  at  each  10  percent 
of  the  blade  span.  The  required  values  of  (AV0/V0t)  and  (AP/AP)  are  calculated 
at  the  appropriate  values  of  x/£  from  the  spline  curves  representing  the 
following  typical  ARL  section  properties: 


X 

£ 

0.0 

0.2 

0.4 

0.8 

1.0 

AV0 

V0t 

0.000 

0.180 

0.427 

0.890 

1.000 

AP 

AP 

0.500 

0.978 

0.171 

1.127 

0.000 

* 

In  PRESKEW,  two  computing  stations  are  used;  however,  only  the  parameters  from 
one  of  these  stations  are  used  to  specify  the  skew  effects  in  SCM-IV. 


f 


■ \ 


i 
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These  section  properties  were  supplied  to  PRESKEW  in  the  form  of  data  statements. 

Once  the  values  of  CU2  and  Cp  are  calculated  from  data  obtained  from 
the  spline  curves  (Equations  13  and  15),  Fy,  AM0 , and  ACPT  are  computed  from 
Equations  (6),  (9),  and  (10),  respectively. 

The  Fy  versus  Y data  are  supplied  to  the  mainline  program  of  SCM-IV,  , 
as  data  statements  and  are  used  in  the  evaluation  of  Equation  (8) . The  AMg 
and  ACPT  data  versus  Y are  supplied  to  the  R0T04  subroutine  of  SCM-IV,  again, 
in  the  form  of  data  statements.  If  a second  skewed,  rotating  blade  row  is 
present  in  the  flow  field,  a separate  PRESKEW  analysis  is  performed  for  this 
blading.  The  resulting  output  parameters  at  the  designated  computing  station 
are  supplied  to  the  mainline  program  and  to  the  STTR4  subroutine  of  SCM-IV  via 
data  statements.  Thus,  the  radial  distribution  of  the  blade  skew  effects  are 
specified  at  one  computing  station  within  the  blade  row.  The  effects  of  the 
skewed  blading  on  the  neighboring  stations  are  affected  by  the  shifting  of 
the  spline  curves  representing  the  streamline  pattern  and  the  accompanying 
redistribution  of  the  meridional  velocity  and  static  pressure  profiles. 

In  the  following  section,  a sample  calculation  based  on  a recent  wind 
tunnel  test  program  employing  a pumpjet  propulsor  with  a skewed  rotor  is 
presented . 

Sample  Calculation  and  Experimental  Comparison 

The  subject  propulsor  is  shown  schematically  in  Figure.  3.  The  required 
input  data  are  presented  in  tabular  form  to  illustrate  the  input  format  to 
PRESKEW. 


Card  No . 1 


XN  = 13.0 
XL AM  = 1.00 

SKEW  = 0.694 


the  number  of  rotor  blades 

the  blockage  effects  were  neglected 

the  constant  blade  skew  angle  (39.7°) 
in  radians 


AD RAT  = 2.092 

Card  No.  2 

XXX(l)  = 1.389 

XXX(2)  = 1.495 


the  average  advance  ratio  used  in  the 
test  program 


the  dimensionless  axial  distance  from 
the  reference  station  to  the  first 
computing  station 

the  dimensionless  axial  distance  from 
the  reference  station  to  the  second 
computing  station 
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Cards  No.  3 

to  6 

X(I) 

Y(I) 

XX(I) 

YY(I) 

1.204 

0.361 

1.279 

0.710 

1.279 

0.710 

1.654 

0.639 

1.575 

0.282 

1.654 

0.639 

1.204 

0.361 

1.575 

0.282 

Cards  No.  7 

to  16 

3i 

C_h 

0.361 

0.000 

0.400 

0.000 

0.440 

0.000 

0.480 

0.000 

0.520 

0.000 

0.550 

0.000 

0.570 

0.000 

0.590 

0.000 

0.610 

0.000 

0.639 

0.000 

Cards  No.  17 

to  26 

Y_ 

CU„ 

E 

E 

0.282 

0.553 

0.300 

0.523 

0.350 

0.467 

0.400 

0.440 

0.425 

0.420 

0.450 

0.400 

0.500 

0.366 

0.550 

0.348 

0.600 

0.350 

0.639 

0.402 

The  end-point  coordinates  of  the 
four  straight  lines  representing 
the  blade  row. 


The  peripheral  velocity  distribution 
at  the  rotor-inlet  plane.  There  was 
no  preswirl  in  the  flow. 


The  experimental  peripheral  velocity 
distribution  at  the  rotor-exit  plane. 
For  initial  blade  specification,  this 
distribution  would  be  dictated  by 
the  design  conditions. 


i 


... 


* 


Cards  No.  27  to  36 


Y 

REF 

cmref 

ANGL(I) 

0.611 

0.500 

- 

0.209 

0.625 

0.590 

- 

0.209 

0.650 

0.665 

- 

0.209 

0.675 

0.707 

- 

0.209 

0.700 

0.7A5 

- 

0.209 

0.725 

0.777 

- 

0.209 

0.750 

0.805 

- 

0.209 

0.765 

0.822 

- 

0.209 

0.785 

0.8A2 

- 

0.209 

0.810 

0.865 

_ 

0.209 
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The  meridional  velocity  profile 
at  the  reference  station.  This 
is  the  same  reference  profile 
that  will  also  be  used  in  the 
SCM-IV  analysis. 


The  PRESKEW  output  data  for  the  XXX(l)  = 1.389  computing  station,  shown  in 
Figure  A,  were  used  in  the  SCM-IV  analysis. 

Discussion  of  Results 


The  results  of  using  the  PRESKEW  output  data  in  a SCM-IV  analysis  of  the 
pump jet  represented  by  Figure  3 are  graphically  presented  in  Figures  A,  5,  and 
6 for  measurements  and  calculations  at  the  rotor  exit  plane.  In  each  of 
these  figures  three  types  of  data  are  presented: 

(1)  a theoretical  analysis  by  SCM-IV  with  shroud  inlet  and 
blade  secondary  flow  losses  included, 

(2)  a theoretical  analysis  by  SCM-IV  with  the  losses  of  (1) 
plus  skew  effects,  and 

(3)  the  results  of  an  experimental  wind  tunnel  investigation 
utilizing  the  subject  propulsor. 

In  general,  the  calculated  data  including  the  effects  of  blade  skew  agree 
more  closely  with  the  experimental  data  than  does  the  SCM-IV  results  employing 
only  the  inlet  and  secondary  flow  losses.  As  predicted  by  Bruce  and  McBride 
(Appendixes  B and  C) . the  total  pressure  coefficient  is  relatively  unaltered 
by  the  effect  of  blade  skew,  whereas,  the  meridional  velocity  exhibits  a 
retardation  in  the  tip  region  with  an  accompanying  acceleration  near  the  root; 
correspondingly,  there  is  an  increase  in  static  pressure  coefficient  in  the  tip 
region  and  a decrease  in  static  pressure  at  the  root.  These  changes  are 
beneficial  in  reducing  cavitation  in  the  rotor  tip  region  and  in  delaying  flow 
separation  at  the  root  section. 
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As  Is  required,  a variable  skew  option  will  be  added  to  PRESKEW;  and 
as  more  is  understood  concerning  the  streamline  curvature  convergence 
problems  additional  computing  stations  will  be  added  within  the  blade  row. 
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CALCULATE  THE  RADIAL 
DISTRIBUTION  OF  THE  DIMEN- 
SIONLESS RADIAL  BODY  FORCE 
AT  THE  COMPUTING  STATION 


USE  SPLINE  CURVE  PRO- 
CEDURES TO  DETERMINE 
GEOMETRIC,  PRESSURE,  AND 
VELOCITY  PARAMETERS  AT 
EACH  TEN  PERCENT  OF  THE 
RADIAL  DISTANCE  ACROSS 
THE  DUCT  AT  THE  COMPUTING 
STATION 


WRITE  THE  OUTPUT  DATA 


RETURN 


As  used  here,  J is  the  computing  station  index. 
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APPENDIX  B 

DISCUSSION  OF  SKEW  EFFECTS  BY  E.  P.  BRUCE 


The  following  two  documents  by  E.  P.  Bruce  contain  a discussion  of  the 
general  radial  equilibrium  equation  and  explain  how  the  equation  can  be 
applied  within  blade  rows.  These  papers  were  originally  written  as  ARL 
Internal  Memoranda;  and  no  attempt  has  been  made  to  match  their  nomenclature 
with  that  used  in  the  body  of  this  report. 


From:  E.  P.  Bruce 


Sub j ec t : Proof  That  A.  L.  Treaster*s  and  L,  H.  Smith’s  Radial 

Equilibrium  Equations  are  Equivalent 

Reference:  Smith,  L.  H. , Jr.,  "The  Radial  Equilibrium  Equation 

of  Turbomachinery,"  A.S.M.E.  Journal  of  Engineering  for 
Power,  pp.  1-12,  January  1966. 

Abstract:  This  memorandum  presents  a short  algebraic  exercise 

which  proves  that  the  radial  equilibrium  equation 
developed  for  use  in  the  Water  Tunnel  SCM  computer 
program  is  identical  to  a simplified  form  of  the 
general  radial  equilibrium  equation  derived  by 
L.  H.  Smith,  Jr. 


Introduction 


In  Reference  1,  the  general  three-dimensional  unsteady  radial 
equilibrium  equation  of  turbomachinery  is  developed.  A rotating,  cylindrical 
coordinate  system  is  used  and  the  resulting  equation  is  then  compared  with  the 
equation  that  would  have  been  obtained  if  a fixed  coordinate  system  had  been 
used.  In  the  development  of  the  Streamline  Curvature  Method  (SCM)  computer 
program  at  the  Water  Tunnel,  A1  Treaster  used  a fixed  coordinate  system.  Thus 
the  two  independently  derived  equations  should  be  equivalent.  Since  the 
equations  are  presented  in  different  forms,  their  equivalence  is  not  immediately 
obvious.  The  purpose  of  this  memorandum  is  to  show  that  the  equations  are 
equivalent. 


Smith* s Equation 


The  equation  of  interest  is  Eq.  (31)  of  Reference  1.  This  is  an  exact 
form  of  the  axisymmetric  radial  equilibrium  equation.  With  Mach  Numbers  set 
equal  to  zero,  this  equation  is 


1 3p 
p 3r 


3 

sec 


C tan  <f> 
z 


(r  tan  <}>) 
3r 


r 
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where 


fluid  mass  density 
fluid  static  pressure 


radial  coordinate 


circumferential  component  of  the  absolute  velocity 
meridional  angle 

axial  component  of  the  absolute  velocity 
radius  of  curvature  of  meridional  pathline 


Treaster's  Equation 

This  equation  has  the  form 


2 2 

. - c C 3C 

— = — + cos  <}>  — - C sin  <p  -r— — 

p 3r  r r m dm 

m 


where  all  quantities  are  as  defined  above  and  Cm  is  the  meridional  component 
of  the  absolute  velocity. 

Proof  of  Equivalence 
2 

The  Cu/r  term  is  common  to  both  equations.  Hence  it  is  necessary  to 
show  that 

C2  3C  , C2  C2  tan  <(>  , _ .. 

, m „ . , m 3 . z , z (r  tan  <f>)  ,,  x 

cos  * — “ Cm  Sin  * 35T  sec  * 7“  + 7 "dr  " (1) 

m m 


Two  relations  are  required.  The  first  is 


C = C /cos  <p 
m z 


and  the  second  can  be  obtained  from  Smith's  Eq.  (81)  which  relates  the  variation 
of  meridional  streamtube  area,  A , to  variations  in  meridional  streamtube,  <p  , 
for  an  axisymmetric  flow  field.  This  equation  is 


1 dA  * d(r  tan  <J>)  . tan  <J) 

A dm  rdr  r_ 


k 
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Along  a streamtube,  the  volume  flow  rate*  Q > is  constant  and  is  given  by 

Q - A C 

m 

From  this 


and  then 


1 3A  C»  3(Q/Cm>  _ 1 9C„ 


A 3m  Q 3m 


C 3m 
m 


3C 


m 


3m 


- C 


m 


cos 


3(r  tan  <{>)  + tan  <{> 
r3r 


m 


Substituting  the  expressions  for  and  3Cm/3m  into  the  left  hand  side  of 
Eq.  (1)  yields  the  following: 


C 3C 

. m - . , ra  _ cos 

cos  <f>  — - Cm  sin  4)  ^ ~ 

m m 


2 A 

cos  <p 


+ sin  4> 


2 A 
COS  (J) 


, 3(r  tan  d>)  , tan  4> 
cos  (0  

in 
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From:  E.  P. 


Bruce 


Subject:  The  Form  of  the  Radial  Equilibrium  Equation  for  Application 

Inside  Skewed  Blade  Rows 

References:  1)  Bruce,  E.  P.,  "Proof  That  A.  L.  Treaster's  and  L.  H.  Smith's 

Radial  Equilibrium  Equations  are  Equivalent,"  ORL  Internal 
Memorandum  File  No.  72-207,  October  12,  1972. 

2)  Smith,  L.  H. , Jr.,  "The  Radial  Equilibrium  Equation  of 
Turbomachinery,"  A.S.M.E.  Journal  of  Engineering  for 
Power,  pp.  1-12,  January  1966. 


Abstract:  A form  of  the  radial  equilibrium  equation  which  can  be  applied 

in  flow  regions  inside  blade  rows  is  presented.  The  new  terms 
of  the  equation  are  discussed  and  their  relation  to  known  flow 
and  geometric  properties  is  defined. 


Introduction 


In  reference  1,  the  radial  equilibrium  equation  developed  for  use  in  the 
SCM  computer  program  is  shown  to  be  identical  to  a simplified  form  of  the 
general  radial  equilibrium  equation  developed  in  Reference  2.  Both  of  these 
equations  are  applicable  to  flow  regions  outside  of  blade  rows.  Additional 
terms  appear  in  the  equation  if  it  is  applied  to  flow  regions  inside  blade 
rows.  The  purpose  of  this  memorandum  is  to  present  and  discuss  the  equation 
containing  these  new  terms. 

The  Radial  Equilibrium  Equation  for  Applications  Within  Blade  Rows 

The  pertinent  equation  is  Equation  (51)  of  Reference  2.  This  equation 
has  the  following  form  when  terms  involving  products  of  Mach  numbers  and 
heat  addition  are  neglected 


i a C 

1 3p  _ u 
p 3r  r 


z 


+ W W 
r z 


3(r  tan  ({>)  + JL  DA 
r3r  X Dz 


+ F 

r 


(1) 


where 

p = fluid  mass  density 

p = fluid  static  pressure 

r ■ radial  coordinate 

Cu  = circumferential  component  of  the  absolute  velocity 
D = differential  change  following  particle 


I 

i 
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z = axial  coordinate 

W = axial  component  of  the  relative  velocity 
z 

Wr  ■ radial  component  of  the  relative  velocity 

$ ■ meridional  angle 

X = effective  area  coefficient 

= radial  component  of  the  body  force  per  unit  mass 


t 


2 sec  <J> 

Setting  W = W tan  <p  , W = C , — 5-  = — ^ (Eq.  37  of  Reference  2)  and 

r z z z * 2 2 

Dz  W r 
z m 

Wm  = Cz  sec  4>  permits  rewriting  Eq.  (1)  as 


I3e- jl  + 

p 3r  r 


3 z 2 
sec  d>  — + C tan  <b 
Y r z 
m 


9(r  tan  ()>)  1.  DX 

r3r  X Dz 


+ F 


Using  the  results  of  Reference  1,  this  equation  may  be  expressed  as 


1 dp  u . . m j a 

— — = — + cos  <p C sin  <p  » — . — 

pdrr  r m T9m  X 


8C  C tan  <{) 

m + _z DA 

9m  X Dz 


(2) 


This  equation,  with  the  omission  of  the  last  two  terms,  is  the  radial 
equilibrium  equation  presently  used  in  the  SCM  computer  program.  In  the 
last  two  terms,  X , DA/Dz  and  Fr  represent  quantities  that  are  not  now 
included  or  derivable  within  the  SCM  program. 

The  effective  area  coefficient,  X , is  given  by 


X 


where 

N = number  of  blades 

0 = circumferential  coordinate  of  blade  suction  surface 

s 

0p  = circumferential  coordinate  of  adjacent  blade  pressure  surface 


j 


Thus  at  a fixed  radius,  a plot  of  X versus  axial  distance  z will  have  the 


where  an  allowance  has  been  made  for  blade  boundary  layer  displacement  thickness, 
i.e.  W 1 at  the  blade  trailing  edge.  This  sketch  shows  the  importance  of 
taking  a number  of  computing  stations  within  the  blade  row.  If,  for  example, 
only  one  station  near  the  center  of  the  blade  row  was  taken,  DA/Dz  ~ 0 and 

2 

Cz  tan  ^ DA 

the  effect  of  the  term ; — — would  be  lost. 

A Dz 

The  radial  component  of  the  body  force  per  unit  mass,  Fr  , is  easily 
evaluated  with  the  aid  of  the  following  sketch: 


The  force,  F'  , on  small  element  of  area  having  dimensions  ds  along  the  blade 
surface  and  dz  axially  is  given  by 


F'  = Ap  ds  dz 
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where  Ap  Is  the  difference  between  the  static  pressure  on  the  pressure  and 
suction  faces,  respectively.  From  the  sketch. 


F*  ■ Ap  ds  dz  sin  £ 


The  radial  force  per  unit  mass  is  given  by 


F = F'/P  dr  r (0  - 6 )dz 

r r v s p' 


Ap  sin  £ 


m I ds 


pr(6s  - 6p)  \dr 


Now  dr/ds  = cos  E so  we  have 
m 


Ap  tan  £ 


F = — 

r pr(6  - 0 ) 
s p 


In  our  designs,  the  quantities  Ap,  tan  and  (0S  - 0p)  will,  in  general,  be 
functions  of  the  radius  r and  of  the  axial  coordinate  z. 

Summary 


A form  of  the  radial  equilibrium  equation  which  can  be  applied  within 
blade  rcws  has  been  presented  and  discussed.  Inclusion  of  this  equation  in 
our  SCM  computer  program  will  require  that  a preliminary  blade  design  be 
completed  which  can  be  used  to  specify  initial  estimates  of  flow  and  geometric 
properties  within  the  blade  passage.  The  importance  of  taking  a number  of 
computing  stations  within  the  blade  passage  has  also  been  noted. 
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APPENDIX  C 

DISCUSSION  OF  SKEW  EFFECTS  BY  M.  W.  McBRIDE 


This  document  which  was  initially  presented  as  an  ARL  Internal  Memorandum 
presents  a discussion  of  the  physical  aspects  of  blade  skew.  The  nomenclature 
is  independent  from  that  used  in  the  body  of  this  report. 


i 


From:  M.  W.  McBride 

Subject:  Geometric  Interpretation  of  the  Effects  of  Blade  Skew 


Reference:  (1)  Bruce,  E.  P.,  "The  Form  of  the  Radial  Equilibrium 

Equation  for  Application  Inside  Skewed  Blade  Rows," 
ORL  Internal  Memorandum  72-211,  October  17,  1972 

Abstract:  The  effect  of  blade  skew  has  been  shown  previously  in 

Ref.  (1).  This  report  intensifies  the  physical  concept 
of  skew  effects  and  arrives  at  a relation  for  the  radial 
pressure  gradient  due  to  skew  in  a blade  row.  This 
relation  is  identical  to  that  determined  in  Ref.  (1). 


Introduction 

The  effect  of  backward  skew  in  a blade  row  is  to  increase  the  static 
pressure  at  the  blade  tip  and  decrease  the  pressure  at  the  root.  This  in 
turn  decreases  the  tip  meridional  velocity  and  increases  the  velocity  at 
the  root.  The  radial  equilibrium  equation  derived  in  Ref.  (1)  clearly 
indicates  this  effect,  but  does  not  give  a clear  physical  indication  of 
the  actual  mechanism  involved.  A simple  geometric  analysis  is  presented 
in  this  report  which  allows  an  understanding  of  skew  while  arriving  at  a 
radial  pressure  gradient  term  identical  to  that  determined  in  Ref.  (1). 

The  Geometric  Analysis  of  Blade  Skew 

A radial  plane  through  a blade  row  is  characterized  by  N evenly  spaced 
blades  with  arbitrary  spanwise  loading  and  thickness  distributions  and 
some  skew  angle  as  a function  of  the  radius.  At  a radius,  R , the  blade 
properties  are: 

a blade  pressure  difference  between  the  pressure  and  suction 

surface  AP; 

a blade  thickness  in  the  tangential  direction  T; 

a skew  angle  £» 

a blade-to-blade  angular  spacing,  o , equal  to  2it/N. 


These  quantities  are  shown  in  the  following  figure.  This  analysis  is  for  the 


■ 
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skew  effect  only  and  neglects  the  pressure  gradient  due  to  rotation  and 
meridional  curvature.  These  terms  are  simply  added  to  determine  the  full 


Assuming  a linear  variation,  the  tangential  pressure  gradient  at  radius  R 
from  blade  suction  surface  to  the  next  blade  pressure  surface  is 


J 

4 


— = AP/(Ra  - T)  . 

ds(R) 


At  radius  R + dr 


—•  = AP/[(R  + dr)  a - T]  . 

ds(R  + dr) 


A value  for  the  pressure  a distance  s from  the  blade  suction  surface 
at  a radius  R,  can  be  determined; 

pl  = p0  + s • Ap/(Ra  - T)  . 

The  pressure  along  a line  eminating  radially  from  that  previously  found, 
but  at  a distance  R + dr  is. 
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k 


?2  «=  PQ  + ( s - dr  tan  C)  AP[(R  + dr)  a - T]  . 

The  difference  between  these  two  pressures  will  give  the  radial  pressure 
change  that  the  blades  are  exerting  on  the  fluid; 

dP  = ?2  - P = (s  - dr  tan  £)  AP/[(R  + dr)  a - T]  - AP/(Ra  - T)  . 

Dropping  the  second  order  terms: 

dP  = - dr  tan  £ AP/(Ra  - T) 


dP  _ AP  tan  g 
drB  " ” Ra  - T * 

The  pressure  gradient  in  the  fluid  is  then  the  negative  of  the  quantity 
dP/drg,  which  is  identical  to  the  term  Fr  defined  in  Ref.  (1).  This  gradient 
is  seen  to  be  uniform  everywhere  between  the  blades  and  is  unaffected  by 
crossing  a blade  surface,  although  the  magnitude  of  the  pressure  will  change. 

- 


j 


I 


w 
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For  blades  whose  skew  angle,  £ , is  invariant  with  radial  distance  the  mean 
thickness  line  must  follow  a logarithmic  spiral  whose  equation  is  given  by 


R 


etan  K 


where 


the  hub  radius,  ft 


= the  peripheral  coordinate,  radians 


FIGURE  1.  The  Basic  Geometry  of  Skewed  Blading 


i 


i 


k.. 


Experimental  Traverse 


.FIGURE  3.  Pumpjet  Geometry-! 


FIGURE  A.  PRESKEW  Data  at  the  XXX(l) =1.389  Computing  Station 


QNV  rldOV  'GWV 


FIGURE  5.  Comparison  of  Experimental  and  Calculated  Data:  Meridional  Velocity  Ratio 
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FIGURE  6.  Comparison  of  Experimental  and  Calculated  Data:  Static  Pressure  Coefficient 
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